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Zero temperature lattice Weinberg - Salam model for the values of the cutoff A ~ 1 
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The lattice Weinberg - Salam model at zero temperature is investigated numerically. We consider 
the model for the following values of the coupling constants: the Weinberg angle 6w ~ 30°, the fine 
structure constant a ^ the Higgs mass Mh ^ 150 GeV. We find that the fluctuational region 
begins at the values of the cutoff A above about 0.8 TeV. In this region the average distance between 
Nambu monopoles is close to their sizes. At A > 1.1 TeV the Nambu monopole currents percolate. 
Within the fluctuational region the considered definitions of the scalar field condensate give values 
■ that differ from the expected one 2Mz/gz- We consider the given results as an indication that the 

nonperturbative effects may be present in the Weinberg - Salam model at the large values of the 
' cutoff. Our numerical results were obtained on the lattices of sizes up to 16'^ x 32. 
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I. INTRODUCTION 
> . 

, Investigation of the phase transitions often requires application of nonperturbative methods. In particular, the 
^ I nonperturbative phenomena are important for the description of the finite temperature Electroweak phase transition 
D , p]-[l4j. At the same time, the phase diagram of the lattice Weinberg-Salam model at zero temperature also contains 
. the phase transition fl5l-[l8|. The phase transition surface separates the Higgs phase from the symmetric phase. On 
both sides of the phase diagram it is necessary to find the way to approach continuum physics within the lattice 
model. It is expected, that the continuum physics arises in some vicinity of this transition (on different sides of the 
^ " transition different continuum models appear). Strictly speaking, this pattern is self - consistent only if the transition 
. is of the second order. In this case close to the phase transition lattice spacing tends to zero and a kind of a continuum 

■ field theory appears. Indeed, zeroth order of the perturbation theory predicts the second order phase transition in the 
] Weinberg - Salam model. However, already on the one - loop level the Coleman - Weinberg effective potential predicts 

the first order phase transition. This discrepancy points out to numerical lattice methods as to the judge. There may 
OO ' take place the 1-st order phase transition or the 2-nd order phase transition. Also another possibility appears: the 
, transition may appear to be a crossover. 

That's why we expect nonperturbative effects to appear in the lattice Weinberg - Salam model close to the transition 
between the two above mentioned phases. Or, in the other words, we expect nonperturbative effects to appear above 
some energy scale, because the increase of the energy scale corresponds to the decrease of the lattice spacing and, 
therefore, is achieved when the phase transition is approached. Basing on trivial dimensional analysis we may expect 
that the mentioned scale can be compared to the Electroweak scale ~ 250 GeV. In fact, some indications were recently 

■ found that this scale might be around 1 TeV (see, for example, [l9l - [25j ). Namely, in the Electroweak theory there 
exist objects that are not described by the first orders of the perturbation theory: Nambu monopoles and the Z - 
strings [2^ . It has been found that there exists the vicinity of the phase transition , where the average distance 
between the Nambu monopoles is compared to their sizes. This region was called in (2^ the fluctuational region 
(FR). Nambu monopoles may be considered as embryos of the unphysical phase within the physical one. Therefore, 
it is natural to suppose that within the FR both phases are mixed and neither the perturbation expansion around 
vacuum with zero scalar field nor the perturbation expansion around vacuum with nonzero scalar field cannot give 
the correct description of the situation. Besides, in [25] it was shown that there exist different ways to define scalar 
field condensate that give identical values out of the FR, and different values within the FR. On the boundary of the 
FR the lattice spacing a remains finite and practically does not depend on the lattice size (for the considered lattices). 
Actually, the value of the ultraviolet cutoff ^ on this boundary for the values of the Higgs boson mass 100, 150, 300 
GeV is around 1 TeV. Thus the hypothesis was suggested that above the energy scale 1 TeV in the Weinberg - Salam 
model nonperturbative effects may become important. It is worth mentioning that these effects, most likely, are 
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related to the expansion in powers of A while the first orders of the perturbation theory for the expansion in powers 
of a are expected to stay at work. In particular, no discrepancy was found between the renormalized fine structure 
constant and its one - loop estimate [25l |. 

In the present paper we extend the research of ^25"] to larger lattices (in [25*1 the lattices of sizes 8'^ x 16 were 
used; here we use lattices 16"^ x 32). In addition wc investigate the properties of Nambu monopoles and Z - strings 
that were out of the scope of the mentioned above papers. Namely, we consider their percolation properties that are 
related to their possible condensation. In the finite temperature theory it was found that the Electroweak transition 
is accompanied with the condensation of the Nambu monopoles and the condensation of the Z - strings [28-30]. Here 
we find that in the zero - temperature model this occurs as well. The "percolation transition" in the Weinberg - 
Salam model at the values of couplings we consider is situated within the FR. And there exists the subregion of the 
FR, were both Z - strings and Nambu monopoles are condensed, while at least one of the considered definitions of 
the scalar field condensate still gives nonzero value. 

The paper is organized as follows. In Section 2 we consider the definition of the lattice regularized Weinberg - 
Salam model and describe the details of the simulation. In Section 3 we discuss the phase diagram of the lattice 
model and the lines of constant physics. In Section 4 we describe how lattice spacing was calculated in our study. 
In Section 5 we calculate the renormalized fine structure constant. In Section 6 we investigate three different scalar 
field condensates. In Section 7 we calculate Z - string and Nambu monopole percolation probabilities. In Section 8 
we discuss the obtained numerical results. Throughout the paper the notations of differential forms on the lattice are 
used (for their definition see, for example, (27|). 



II. THE LATTICE REGULARIZED WEINBERG - SALAM MODEL 



We consider the model without fermions. Its partition function has the form: 



Z = J DHDrexp{-A{r,H)) (1) 
Here ^(F, H) is the action for the scalar doublet H and the gauge field T ^ U <^ e''^ e SU{2) » [/(I): 
A{r,H) = P ((l-5Tr[/p) + -^(l-cos0p)) + 

plaqucttcs 

7 ■ 



+ ^ E - U.ye^'-'''Hy\^ + Y.m.?{l - 2A - 47) + X^'), (2) 

xy X 



On the tree level we have: 



,2^ 

'8A 

7 



7 

mz = V 



7c = (3) 

Here we have introduced: vacuum expectation value v of {H^l, the lattice Higgs boson mass rriH ~ Mhcl, the lattice 
Z - boson mass mz = Mza, and the critical value 7c. 
After fixing Unitary gauge 

i/= J) (4) 

where H is the scalar doublet, the Z2 gauge ambiguity remains: — > (— l)"^/!^;, Z ^ [Z + 7rdri]mod 27r. Here the Z 
- boson field is defined as 

Z = -Arg [Une''] (5) 
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FIG. 1: The phase diagram of the model in the (7, A)-plane at /3 = 12. The dashed line is the tree - level estimate for the line 

of constant physics i:^ = s^'ST^ ~ const) correspondent to bare = 150 Gev. The continuous line is the line of the phase 

transition between the physical Higgs phase and the unphysical symmetric phase (statistical errors for the values of 7 at each 
A on this line are about 0.005). The dotted line is the line A = 0.0025. Along this line the physical quantities are calculated in 
the present research. 



The tree level approximation gives for the infrared effective constraint potential [25| after any Z2 gauge is fixed: 

V'-'-i^) = NiX{^^-v^f (6) 

Here A'4 is the lattice volume. 

In numerical simulations we use Metropolis algorithm. The model is simulated in Unitary gauge with the signs of 
h unfixed. After each 150 Metropolis sweeps the Z (or DZ) - version of Unitary gauge is fixed (for the definition of 
these gauges see the next sections). As a starting point of our simulations on the lattice 16'^ x 32 we use configurations 
obtained on the lattice 8'^ x 16 during the preparation of [25j. 16 identical configurations are merged together forming 
the starting 16'^ x 32 configuration. Then, about 60000 Metropolis sweeps are made before the measurement of 
observables begins (this has required about 600 ours CPU time). During this preliminary run the 16 mentioned above 
parts of the lattice become decorrelated which signalizes that the thcrmalization is achieved. 



III. PHASE DIAGRAM AND LINES OF CONSTANT PHYSICS 



The lattice model defined by Eq. ^ has the four - dimensional {l3,j,X,9w) phase diagram. On this phase 
diagram phase transition surface is three - dimensional. The lines of constant physics on the tree level are the lines 

i:^ = g^^T- = const; /? — — const; 9w — const). We suppose that in the small vicinity of the transition the 
deviation of the lines of constant physics from the tree level estimate may be significant. However, qualitatively their 
behavior is the same. Namely, the cutoff is increased along the line of constant physics when 7 is decreased and 
the maximal value of the cutoff is achieved at the transition point. Nambu monopole density in lattice units is also 
increased when the ultraviolet cutoff is increased. 

In our lattice study we fix bare 9w = ""/B,/? — 12, A = 0.0025. Therefore, strictly speaking we investigate the 
system along the line on the phase diagram that differs from the line of constant physics. This is illustrated by Fig. 
[1] where the projection of the phase diagram to the plane (/3 — 12, 9w = ""/S) is drawn. This diagram is obtained, 
mainly, using the lattice 8'^ x 16. Some regions (A = 0.009, 0.0025, 0.001), however, were checked using larger lattices 
(see, for example, [l9l. [20l. [22| - [23 | ) . According to our data there is no dependence of the diagram on the lattice size. 
The physical Higgs phase is situated right to the solid transition line. The position of this line was determined using 
various methods 19, 20, 22" 24]. However, the uncertainty is still present in the final determination of the phase 
transition points. This uncertainty on this Figure is within the error bars of 7 = 7c ± 0.005. The details of this 
uncertainty at A = 0.0025 are discussed in the present paper (see below, this Section). The dashed line represents 
the tree - level estimate for the line of constant physics. The dotted line is the line A = 0.0025. Along this line the 
physical quantities are calculated that are reported in the present paper. One can see, that already on the tree level 
this straight line deviates from the line of constant physics. However, Figure [T] demonstrates also that within the 
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FIG. 2: The link part of the action as a function of 7 at A = 0.0025 , /? = 12 on the lattice 16^ x 32. 



interval 7 G [0.255, 0.27] (where the physical quantities are measured in the present research) the deviation of the 
tree level estimate for the line of constant physics from the straight line A = 0.0025 is not crucial. In fact, on the 
tree level along this straight line the fine structure constant does not vary. The renormalized fine structure constant 
is also almost not changed (see discussion below in Section |V]). As for the Higgs boson mass, its value on the tree 
level varies between 154 GeV at 7 = 0.255 and 145 GeV at 7 = 0.27. The variation of the renormalized Higgs mass 
along the line A — 0.0025 is discussed in Section HVl Our data (with large statistical errors, though) also demonstrate 
that the Higgs mass does not deviate significantly from bare value ~ 150 GeV at A = 0.0025,7 € [0.2585,0.27]. It is 
worth mentioning that we did not investigate the renormalized W - boson mass in the present research. Therefore we 
do not represent here any data on the renormalized Weinberg angle. However, we also expect that it does not vary 
sufficiently at the considered values of 7. 

In the present paper we deal with the line A — 0.0025, 9 = 7r/6, (3 — 12 for 7 S [0.2585, 0.27] as with an approximation 
of the line of constant physics because along this line the main physical quantities (Higgs mass, fine structure constant) 
do not vary essentially. The lowest value 7 = 0.2585 in the mentioned above interval is chosen because we expect that 
for the description of the model at 7 < 0.2585 larger lattice sizes are necessary. This is related to the fact that at 
7 < 0.2585 the value of the cutoff is larger than 1.5 TeV and increases very fast (see discussion below, in Section [TVl 
Eq. pip ). At the same time already for the values of the cutoff A ^ 10 TeV we need lattices of linear size >> 10 (see 
discussion in Section rVIII[) . We expect that the line of constant physics at A >> 1 TeV deviates essentially from the 
straight line investigated in the present paper. Not only the Higgs mass may deviate from its bare value but also the 
renormalized fine structure constant and, probably, the renormalized Weinberg angle. However, for the investigation 
of such large values of A extremely large lattices are needed and such a research is out of the scope of the present 
paper. 

Below in the present paper we always deal with the investigation of the lattice model at fixed 9w = tt/G,/? = 
12, A — 0.0025. In Fig. [2] the data of the the link part of the action j^J^xy^xUxye^^^^Hy are represented. The 
dependence of the link part of the action on 7 indicates that the phase transition (or, a crossover) can be localized at 
7 ~ 7c = 0.25775 ± 0.00025. (See also [2^, where the same conclusion was made basing on the data obtained on the 
lattice 8^ X 16.) As in we exclude the first order phase transition because we do not observe any sign of a two - 
state signal. In the next sections we shall demonstrate that the infrared UZ potential for the scalar field also points 
out to 7^ as to the transition point. 

In addition in Fig. |3]we represent the fluctuation of the zero momentum scalar field x — [SHp^o]'^ =< [Hp^o]'^ > 
— < Hp^o >^ in lattice units as a function of 7. Here i/p=o = From this plot we obtain the critical 

value 7" — 0.258 ± 0.0005. Moreover, the given plot indicates that we deal with the second order phase transition 
localized at 7". Indeed, the values of the fluctuation SHp^o on the lattice S'^ x 16 are about = 4 times larger 
than on the lattice 16^ x 32 for 7 > 7". This means that for these values of 7 the physical size of the lattice 8^ x 16 
is twice smaller than that of the lattice 16'^ x 32, as it should when the correlation length is smaller than the lattice 
size. However, at 7" the fluctuations calculated using both lattices almost coincide with each other. This means that 
at 7 = 7" the physical sizes of both lattices are the same, that may happen only if the correlation length is much 
larger than the lattice size. It is worth mentioning, that the difference between 7^ and 7" is 0.00025 that is within the 
statistical errors of both quantities. We feel this appropriate to refer to 7^ as to the possible phase transition point. 
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FIG. 3: Susceptibility x (in lattice units) as a function of 7 at /3 = 12, A = 0.0025. Circles correspond to lattice 8^ x 16 while 
squares correspond to lattice 16^ x 32. 




0.255 0.260 0.265 0.270 0.275 

GAMMA 



FIG. 4: Z - boson mass in lattice units as a function of 7 at A = 0.0025 , /3 = 12. Crosses correspond to lattice 8'^ x 16. Circles 
correspond to lattice 16^ x 32. 



It is worth mentioning that Fig. [3] allows to estimate the fluctuation of the scalar field within the fluctuational 
region. On its boundary, at 7 = 0.2625 (see below) we have x ^ 0.001 on the lattice 8"^ x 16. Also we know that 
at this value of 7 the correlation length for the scalar field is about two lattice spacings. Therefore the given lattice 
contains 4^^ x 8 = 512 cubes of the linear size equal to the correlation length. This means that the fluctuation of the 
scalar field within such a cube is 5\H\ ^ * x ~ 0-7 that is to be compared with the average value {\H\) ^ 2. 

The formal requirement for the perturbation theory to be applied is 2 ~ {\H\) >> S\H\ ~ 0.7. It seems that this 
inequality is not satisfied. For 7 < 0.2625 the situation is even worse. For example, at 7 = 0.26 we have {\H\) ^ 2 
while S\H\ - 0.9. 

IV. Z - BOSON MASS, LATTICE SPACING, AND HIGGS BOSON MASS 

For the calculation of the Z - boson mass we use the following definition of the Z boson field: 

Z.,y ^Z^ = - [Arg($+C/,,e'''-ci>,)]. (7) 

Actually this definition of the Z - boson field coincides with the previous one ([5| taken in the version of Unitary gauge 
(III) with nonnegative h. 
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In order to evaluate the mass of the Z-boson we use the correlator 

1 

iV6 



x,y fi 



Here the summation y is over the three "space" components of the four - vectors x and y while a:oj 2/o denote their 
"time" components. N is the lattice length in "space" direction. L is the lattice length in the "time" direction. The 
full lattice 4 - volume is iV4 = x L. 

In order to evaluate the Higgs boson mass we use the correlator: 

^^{(|i/,||i/^|)-(|i/,|)2}^e-™-l"l +e-™-(^-|"l) (9) 

We can roughly evaluate the dependence of the lattice Z - boson mass on the lattice size as follows. In finite 
temperature theory gauge boson thermal masses appear of the order of nig = Mga ^ gTa ~ where T is the 
temperature while Nt is the lattice size in imaginary time direction. Analogy to the finite temperature theory allows 
us to evaluate the finite volume contribution to the Z - boson mass as Amz ~ where N is the linear lattice size. 
For a 1/150, TV ~ 8 we have Amz ^ 0.08 while at iV 16 we expect Amz ^ 0.04. 

In Fig. 21 we represent our data on the Z - boson mass. Our numerical results confirm the results of ^24^. Nonzero 
values of Z - boson mass are obtained at 7 > 0.2585. At the same time for 7 < 0.258 we observe large statistical 
errors for the ZZ correlator. Therefore, in this region of the phase diagram the Z - boson mass cannot be calculated 
and we suppose it vanishes somewhere between 7 = 0.25 and 7 = 7c- 

Taking into account expression ([3]) for the Z - boson mass we use the following fit (jc is changed to 7^ = 0.2575) to 
the data of Fig. [J] : 



y pAcos^'t'iy 

This fit is represented on the plot by the dashed line. It is worth mentioning that 7^ is within the error bars of 7^. 
We find that ([10]) with this value substituted instead of 7c approximates the data better than with 7c = 0.25775 
or 7c = 0.258. This does not mean, that Fig. |4] points out to 7c = 0.2575 as to the transition point instead of 
7c = 0.25775. Instead, this means that there exist also other contributions to the dependence of the Z - boson mass 
on 7 in addition to the tree level estimate with the real critical value of 7 substituted instead of its naive estimate. 

Using the value of lattice Z - boson mass mz we can evaluate the ultraviolet cutoff A = ^ as a function of 7 via the 
relation mz = a x 91 GeV. Further we shall use fit ([TOl) in order to represent our results as a function of A. Namely, 
we use the following representation for A: 



/ /3Acos26'w c,^n\T 
^ ^ "V 27(7 - 0.2575) '^'^^^^^ (^^^ 

In particular, in Fig. [5] we represent the Higgs boson mass calculated on the lattice 8"^ x 16 in physical units (in 
GeV) as a function of the cutoff. Unfortunately, we do not have enough statistics to calculate this mass on the larger 
lattices. We observe that the calculated values of the mass are close to the expected value ~ 150 GeV. The deviation 
is within the statistical errors. 

Our estimate for the ultraviolet cutoff at 7c is ~ 1 TeV. The value of Z - boson mass in lattice units at this point 
is about ~ 0.2. The above mentioned estimate of the finite volume effect is ^ 0.04. Therefore we expect the values 
of the ultraviolet cutoff reported here cannot differ from that of obtained on an ideal infinite lattice by more than 20 
per cent. Thus we give here the following estimate for the cutoff A at 7 = 7c: Ac = 1 ± 0.02 TeV. 

V. RENORMALIZED FINE STRUCTURE CONSTANT 

In order to calculate the renormalized fine structure constant a^j = e^/47r (where e is the electric charge) we use the 
correlator of Polyakov lines for the right-handed external leptons. These lines are placed along the selected direction 
(called below the imaginary time direction). The space - like distance between the lines is denoted by i?. 

C(|x - y\) = (ReH(e^'^<*-')<*-'+i) H(e"^'^<«'')<S''+i'). (12) 
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FIG. 5: Higgs - boson mass in physical units as a function of the cutoff at A = 0.0025 , /3 = 12 on the lattice 8^ x 16. 



The potential is extracted from this correlator as follows 

V{R) = -hogCiR) (13) 

Here L is the size of the lattice in imaginary time direction. 

Due to exchange by virtual photons and Z - bosons one would expect the appearance of the Coulomb and Yukawa 
interactions: 

V(r) = -aB.[llo{r) + ^Unizir)] + const, 

""W - -J^ 2^ + sin2p2/2 + + 

Here N is the lattice size, pi = h — 0, L ~ 1. 

We substitute to ^(T^ the fit to the Z - boson mass represented in (ITU)) . The results are presented in Fig. [Bjand are 
to be compared with the tree level estimate for the fine structure constant a'"-* ^ and the 1 - loop approximation 
(when we assume bare value of a to live at the scale ~ 1 TeV while the renormalized value lives at the Electroweak 
scale Mz): a(i)(Mz/l TeV) - 

The values of the renormalized fine structure constant calculated on the lattice 16'^ x 32 are close to the values 
calculated on the lattice 8'^ x 16 represented in [l^. We observe that the renormalized fine structure constant calculated 
in the mentioned above way is rather close to the one - loop estimate (when the cutoff A in ^'^^(M^/A) is around 1 
TeV. It is worth mentioning, that the present data on the renormalized fine structure constant (and the data of |25|) 
differ from that of reported in ■ In Iz^l we used the potential extracted from the Wilson loops and approximated 
it by the simple 1/R fit. Moreover, in [2J] the exchange by virtual Z - bosons was neglected. Therefore the values 
represented in depend strongly on the lattice size and deviate essentially from the one - loop estimate near to the 
phase transition point. 

VI. SCALAR FIELD CONDENSATE 

In [2^ three different effective constraint potentials were introduced. All them are defined in Unitary gauge 
h 

^ 

exp(-F"-"(0)) = (<5(0-/i,)) (15) 
Here real scalar field is defined on the lattice points x. 



H ~ [ ^ ] with real h. The first one is the ultraviolet potential 
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FIG. 6: The inverse renormalized fine structure constant l/a{Mz /A) as a function of the cutoff A at A = 0.0025 , /3 = 12 on 
the lattice 16'^ x 32. The dashed line is the one - loop estimate for A = 1 TeV. 
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FIG. 7: The scalar field condensate (in GeV) as a function of the cutofl for A = 0.0025 , /3 = 12. Circles correspond to the 
UZ potential, lattice 16^ x 32. Empty squares correspond to the UZ potential, lattice 8'^ x 16. Crosses correspond to the UDZ 
potential, lattice 16^ x 32. Stars correspond to the UDZ potential, lattice 8^ x 16. Triangles correspond to the ultraviolet 
potential, lattice 16^ x 32. Diamonds correspond to the ultraviolet potential, lattice 8"^ x 16. 



In order to define the infrared potential 

X 

(where is the number of lattice points) it is necessary to fix the ambiguity 

hx^ {-ly^h.^,, [Z + 7rdn]mod27r, (17) 

where the Z - boson field is defined in ([5]). 
The first way is minimization of 

E(l-cosZ) min (18) 

links 

with respect to the mentioned Z2 transformations. In 25] this gauge was called the Z - version of Unitary gauge and 
the corresponding effective potential (fT6| is called UZ potential. 
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FIG. 8: The scalar field condensate extracted from the UZ potential in lattice units as a function of 7 for A = 0.0025 , /3 = 12 
on the lattice 16^ x 32. 



The second way to define the Unitary gauge with ^ R is to minimize the divergence of Z with respect to the 
remaining Z2 transformations: 



2 



min (19) 



This gauge is cahed the DZ - version of Unitary gauge and the corresponding effective potential (IT6l) is called UDZ 
potential. 

The three mentioned above effective potentials give three different definitions of the scalar filed condensate. (The 
condensate v is defined as the value of (f>, at which the potential V{(j)) has its minimum). 

In Fig. [7] we represent these three condensates as functions of the cutoff A — ^. We consider the condensates in 

physical units, i.e. we multiply the values expressed in lattice units by — . So, we define w^'*!"* — 

The scalar field condensate (in physical units) has to be renormalized. The renormahzed condensate is usually 
defined as f o''^* ~ I/2 v^^"^^ , where the wave function renormalization constant enters the following approximation 

for the scalar field propagator: 

{[Hp]+H„)-\{H)\'' = ZH fi^ 5 , (20) 

" " 7V4(4E,sinV2 + 4sh2ii^)' ^ ' 

Here ffp = ^ e'f^iJ,, and = N^L. 

Due to the renormalizability the Z - boson mass is related to the scalar field condensate as follows: Mz = 
gzv^^^ /2 = —Yj2gzv'^^'''^^/2, where gz — sfi^2e" renormalized coupling constant. In the perturbation thc- 

ory the deviation of Zh from unity is proportional to the factor ~ alog-p- ^ 0.02 (for A ^ 1 TeV). Therefore, for 
A 1 TeV we expect Zh ^ I- 

That's why the perturbation theory prompts that both the (nonrenormalized) scalar field condensate (represented 
in Fig. [7]) and the renormalized one must be close to the value v^^^^ = 2Mz/gz = sin26'vy Mz /V 47ra ~ 273 GeV 
(this value differs from the conventional value 246 GeV due to the difference in a). We observe that all considered 
condensates approach this value when the cutoff is decreased. However, an essential deviation from this value appears 
at the values of the cutoff ^ 1 TeV. Looking at this plot we also conclude that the condensates extracted from the 
UDZ potential and from the ultraviolet potential represent close quantities 33]. At the same time the UZ potential 
gives different value of the scalar field condensate. At the present moment we do not understand what is the reason 
for such a behavior. It is worth mentioning that the condensate extracted from the UZ potential vanishes at 7^ (see 
FiglS]). We must remember, that the cutoff was not calculated at this point but Fig. [3] indicates that at this point we 
may have A^ ^ 00. (Let us remind also that Figl2] points out to 7^ as to the phase transition point.) The other two 
condensates vanish at Ac ^ 1 TeV i.e. at 7c — 0.26075 ± 0.00005, where the value of the cutoff has been calculated 
explicitly. We shall see in the next section that close to this point of the phase diagram both the Z - string and the 
Nambu monopoles begin to percolate. 
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FIG. 9: The percolation probability for the Z - string (circles) and for the Nambu monopoles (crosses) as a function of the 
cutoff A = f (in GeV) at A = 0.0025 , ^ = 12 on the lattice 16^ x 32. 



The observed behavior of the scalar field condensates calculated in UZ and UDZ - versions of Unitary gauge 
means that the wave function renormalization constant for the scalar field (defined in these gauges) differs from the 
perturbation theory prediction at A > 800 GeV. (We may calculate this constant as Zh — ( ^^p^°^ )^, where yP'^^'^ is 
drawn in Fig 171) 



VII. Z - STRINGS AND NAMBU MONOPOLES 



In this section we use definition ^ of the Z - boson field. The classical solution corresponding to a Z-string should 
be formed around the 2-dimensional topological defect which is represented by the integer-valued field defined on the 
dual lattice S = ^ ([dZ,„od27r ~ dZ). Therefore, S can be treated as the worldsheet of a quantum Z-string 28]. Then, 
the worldlines of quantum Nambu monopoles appear as the boundary of the .Z-string worldsheet: jz = STi. 

The percolation probability of Nambu monopoles is defined as follows. First let us denote the probability that two 
points X, y are connected by the monopole cluster by p{x,y). We may identify it with the following quantity 

(*+*,)= p(a:,y), (21) 

where operator ^1'+ creates the monopole - antimonopole pair at the point x. This identification allows us to calculate 
the lightest monopolium mass ttim, i-e. the mass of the quantum state consisted of the monopole - antimonopole pair 
connected by the Z - string: 

ni(ko-2/o|) = ^ J2 (^^^i') 

xi,X2,xa,yi,y2,y3 

^ y4 (e-™Afl^o-aol ^ g-mM(i-|a;o-yol)') (22) 

Here the lattice size is x L, and it is implied that the mass is calculated in the region of the phase diagram, where 
the condensate of ^> vanishes. In order to calculate this condensate we consider the following quantity: 

^ \m'+m^-y\), (23) 

where x — {xq,xi,X2), and ^3{r) — >■ (r — >■ oo). Thus the condensate is defined through the percolation probability: 
Cmon ~ = limr->.oon3(r). It is worth mentioning that the two quantities ni(r) and Tlsir) give different values 

at r cx) if there is the massless scalar excitation in the spectrum, otherwise ni(r) — ^sir) {r oo). 

In a similar way we are also able to calculate the mass of the lightest excitation of the Z - string. First, the 
probability of the two points to be connected by the Z - string cluster is defined. Next, this probability is related to 
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FIG. 10: The monopolium mass (tiangles) and the hghtest string excitation mass (squares) as a function of the cutoff calculated 
at A = 0.0025 , /3 = 12 on the lattice 8^ x 16 (empty symbols) and 16^ x 32 (dark symbols). 
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FIG. 11: The Nambu monopole density (in lattice units) as a function of the cutoff at A = 0.0025 , /3 = 12 on the lattice 
16^ X 32. 



the two - point correlator of the operators that create the Z - string excitations. After that the mass of the lightest 
excitation is extracted and the percolation probability for the Z - strings is defined. The percolation probabilities for 
the Nambu monopoles and the Z - strings are represented in Fig. [9l We observe that both monopole currents and Z 
- string worldsheets percolate at A > 1.1 TeV. 

The monopolium mass and the lightest Z-string excitation mass calculated are represented in Fig. [TU] as functions 
of the cutoff. One can see, that both these masses decrease when the cutoff is increased. At A > 1.1 TeV we do not 
calculate the mentioned masses because the condensates of the monopolium and of the Z string excitations appear. 

The code for the calculation of the percolation probability was written especially for the investigation reported in 
this paper. It has been tested in several ways. In particular, results of 29] and ,3Qi] were reproduced. 

According to the classical picture the Nambu monopole size is of the order of Mjj^. Therefore, for example, for 
~ 250 GeV and Mh ~ 150 GeV the expected size of the monopole is about two lattice spacings. The monopole 
density around 0.015 means that among about 16 sites there exists 1 site that is occupied by the monopole. Average 
distance between the two monopoles is, therefore, about 2 lattice spacings that is the monopole size. In Fig. [Til (where 
the Nambu monopole density is represented as a function of the cutoff) we observe, that at A > Ac2 ~ 0.8 TeV the 
Nambu monopole density is larg er than 0.015, i.e. the average distance between monopoles is less then the classical 
monopole size. According to 24] this means that at Ac2 we enter the fluctuational region. It is worth mentioning that 
within this region the notion of quantum Nambu monopole differs from the notion of the classical Nambu monopole 
considered in [26|. In particular, the size of the quantum object may be sufficiently less than that of the classical 
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Nambu monopole. 

The original estimate of the Nambu monopole mass was given by Nambu [2^ : 

Mat- — sin5 6lvi/\/— 246GeV- 900GeV (24) 
3e V Mw 

Then, according to [26| the classical energy of monopole - antimonopole pair is 

E = 2Mm~^^, (25) 

where Q — ^ is the monopole charge while / is the average distance between the monopoles in the monopolium. 

We can use this formula to estimate roughly the dependence of the monopolium mass on the cutoff. At small values of A 
the contribution of the l/l term can be neglected (as the monopole density is negligible) and the lightest monopolium 
mass is expected to be close to the value ~ 1.8 TeV. However, when the cutoff is increased, the average distance 
between Nambu monopoles is decreased. Therefore, I in is decreased as well. As a result the monopolium mass 
is decreased. Indeed, we observe this kind of behavior in Fig. [TUl 

It is worth mentioning that the monopolium and the Nambu monopole itself are unstable as classical objects 
[2^. However, when the cutoff is increased, the whole picture of Nambu monopoles and monopolium is changed. 
The operator (defined in (|2ip ) that creates the monopole - antimonopole pair creates actually the quasiparticle, 
the properties of which may differ essentially from the properties of the classical monopolium. In particular, these 
quasiparticles are condensed at A > 1.1 TeV. 



VIII. DISCUSSION 



Let us try to estimate the conditions under which the perturbation theory can be applied to this model. We make 
this estimate in the spirit of volume 5, paragraph 146 of [sH ] (where the similar considerations were used in order to 
estimate the width of the fluctuational region in the finite temperature Ginzburg - Landau model). We are going to 
compare the vacuum average of H with the fluctuations of H within the 4- volume N4, the linear size of which is 
equal to the correlation length of H. This correlation length in lattice units is equal to l/rrijj. The fluctuations are 

obtained from JSJ and are of the order ofSH ~ — , ~ Therefore, at u >> '"^ the perturbation theory can 

' •uv'8-'V4A vVSX vVSX _ 

be applied while otherwise it might not be applied. Finally, we obtain A << 7^/8 or ^^^vAna << 1. For mn ~ 150 
GeV and a ~ 1/150 this expression reads 0.54 << 1. This estimate is indeed conflrmed by numerical results (see the 
end of Section 2). Thus already on this level there may appear some doubts about the validity of the perturbation 
expansion within the given model. 

Above we have reported the results of our numerical investigation of lattice Weinberg - Salam model at /3 = 12, 
A = 0.0025, 6w = 30°. For these values of couplings the bare Higgs boson mass is close to 150 GeV near to the 
transition between the Higgs phase and the symmetric phase. Numerical simulations were performed on the lattices 
of sizes up to 16'^ x 32. 

Our data draw the following pattern of the phase transition. 

1. When the cutoff is increased (7 is decreased) Z vortices become more and more dense. Somewhere around 
A ^ 0.8 TeV (7c2 ~ 0.2625) the transition to the fluctuational region occurs [IJ]. In this region Z - vortices and the 
Nambu monopoles dominate. The average distance between Nambu monopoles becomes compared to their sizes. 

2. At the value of A around Ac ~ 1 TeV (7 around 7c ~ 0.26075) the scalar field condensates calculated using 
the UDZ effective potential and the ultraviolet effective potential vanish. At A > 1.1 Tev the Nambu monopoles and 
the Z - strings begin to percolate. This means, in particular, that the operator that (naively) creates the so - called 
monopolium state actually creates the quasiparticles that are condensed. 

3. At the value of 7 around 7c ~ 0.25775 the scalar field condensate calculated using the UZ effective potential 
vanishes. Also somewhere close to 7c the derivative of the link part of the action has the step - like discontinuity. We 
cannot calculate the ultraviolet cutoff at 7c due to large statistical errors. At the present moment we do not exclude 
that it tends to infinity at this point on the ideal infinite lattice. 

The technical question about the order of the phase transition remains. There still exist two possibilities: either 
we deal with the second order phase transition (localized at 7c) or with the crossover. The first possibility is realized 
if at 7c all lattice masses vanish or the correlation lengths are infinite at this point. The behavior of the scalar field 
fluctuation indicates that this may indeed be true. However, accurate investigation of the lattice masses is still to 
be performed in the vicinity of 7c. Let us now estimate the lattice size needed for such an investigation. Suppose, 
we need to investigate the region of the phase diagram with A = ^ 1 TeV. Then, the lattice size has to be much 
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larger than the correlation length: L » = q^v a ^ ^- Therefore, the lattices of sizes 16^ x 32 seem to us 
large enough to investigate the model at 7 > 0.2585. If, however, we are going to investigate the region of the phase 
diagram with A ^ 10 TeV, we need to have lattices of sizes L >> 10. For this purpose lattices used in the present 
research are not large enough. That's why if the second order phase transition is indeed present at 7^, in a small 
vicinity of this point, where A >> 1 TeV (most likely, this vicinity is situated within the interval [0.2575,0.2585)), 
the numerical lattice methods that use lattices of sizes up to IG'^ x 32 cannot be applied for the calculation of lattice 
masses. 

The scalar field condensates calculated in three different ways in our study differ from each other and from the 
expected value ~ 273 GeV for the values of A > 800 GeV. However, all them have tendencies to approach this value 
when A is decreased. This means that the wave function renormalization constant for the scalar field differs from 
its perturbative estimate at A > 800 GeV. The percolation probability for the Nambu monopoles and for the Z - 
strings differ from zero at A > 1 TeV. We consider this behavior as a manifestation of the nonperturbative effects 
present in the given model at large enough energy scales. It is worth mentioning that the point of view that the 
nonperturbative effects may become important in the Higgs sector of the Standard Model is not new. In particular, in 
[s^] it was argued that the wave function renormalization constant for the scalar field contains large nonperturbative 
contribution (at zero temperature). This conclusion of [3^ is in accordance with our results represented here. 

The situation seems to us similar to the phase transition in the second order superconductors at finite tempera- 
ture. Namely, there exists the fluctuational region around the critical temperature {Tc — AT;Tc + AT), where the 
perturbation theory cannot be applied and the nonperturbative effects are present jsij. In the lattice Weinberg - 
Salam model (at zero temperature) there also exists such a region around the phase transition at 7^. This region is 
localized within the interval [7^; 702) that corresponds to values of the cutoff A > 0.8 TeV. As it was explained above 
we observe indications that within this region nonperturbative effects become important. 
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